LINEARIZING NON-LINEAR INVERSE PROBLEMS 
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Abstract. We propose an abstract approach to prove local uniqueness and conditional Holder 
stability to non-linear inverse problems by linearization. The main condition is that, in addition 
to the injectivity of the linearization A, we need a stability estimate for A as well. That condition 
is satisfied in particular, if A* A is an elliptic pseudo-differential operator. We apply this scheme 
to show uniqueness and Holder stability for the inverse backscattering problem for the acoustic 
equation near a constant sound speed. 



1. Introduction 

Non-linear inverse problems are often linearized, and reduced to the problem of the injectivity 
of their linearization. If that linearization is an injective map with a closed range, then this 
implies local uniqueness and Lipshitz stability, see Theorem [1] below. Most inverse problems in 
Mathematical Physics however are ill-posed and do not fall in this category. The mere fact that 
the linearization is injective, if it is, is not enough for local injectivity of the original problem when 
the closed range condition fails. The closed range condition is also equivalent to a linear stability 
estimate, see (fTTj) below. The main purpose of this work is to propose a systematic approach to treat 
inverse problems for which the closed range condition may still hold if we replace the original spaces 
with new ones satisfying interpolation estimates, see (jlip . A typical example is to replace C k or H s 
by C k , respectively H s , with different k' , s'. A sufficient condition for the linear stability estimate, 
besides the injectivity of the linearization A, is A* A to be an elliptic pseudo-differential operator 
(^DO). The later is a standard consequence of the theory of elliptic ^DOs and elliptic 'I'DOs. Our 
main result in this direction is Theorem [2] that states, roughly speaking, that linearization plus an 
appropriate stability estimate of the linearized map imply local injectivity and a conditional Holder 
stability for the non-linear map. 

We recently used this approach for studying boundary rigidity /lens rigidity questions for compact 
Riemannian manifolds with boundary 112 } [T3]. Some of these ideas have been used before, on a 
case by case basis, for proving local uniqueness with or without a Holder stability, see for example 
[3j [91 [10]. Our goal is to systematize this approach, and in particular to understand what makes 
(at least a large enough class of) ill-posed inverse problems ill posed, and in what cases one can 
get local uniqueness by linearization. In particular, we show why many severely ill-posed inverse 
problems cannot be treated this way — the linearization is unstable in any pair of Sobolev spaces. 
The most famous example of this type is Calderon's problem [2l I16|. [T9] . 

We give an application to the inverse backscattering problem for the acoustic equation. In this 
problem, we are trying to determine the sound speed or, equivalently, the index of refraction, 
of a medium by measuring the scattering amplitude with the direction of incidence of a plane 
wave opposite to the direction of the reflected wave. In other words, we are measuring the echos 
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produced by plane waves. This problem arises for instance in ultrasound tomography. We prove 
injectivity near constant sound speeds and we give a Holder type conditional stability estimate, see 
Theorem HJ Conditional Lipshitz stability estimates in different norms were given in [20J based on 
the local uniqueness proof in [10J. 

In section [2] we prove the main result — Theorem [2j In section [31 we give sufficient conditions 
for the stability of the linearized inverse problem required by Theorem [2j Those conditions are 
based on standard elliptic ^DO theory and Fredholm theory. Finally, in section [H we give the 
application to the inverse backscattering problem. 

2. An abstract inverse problem 

We describe our abstract inverse problem now. Let A : B\ — > B2 be a continuous non-linear map 
between two Banach spaces. Consider the "inverse problem": 

(1) Given h <E Ran (.A), find / so that A(f) = h. 
Here, we study only the (local) uniqueness and stability questions for ([I]). 

Definition 1. We say that there is a weak local uniqueness for {7]) near fo € B\, if there exists a 
neighborhood of fo so that for any f in that neighborhood that solves 

(2) A(f) = A(f ), 
we have f = fo- 

Definition 2. We say that there is a strong local uniqueness for |7J) near fo € B2, if there exists 
a neighborhood of fo so that for any f\, $2 in that neighborhood with 

(3) A(h) = A(f 2 ), 
we have f\ = fa. 

Besides uniqueness, we are often interested in stability estimates. By a classical argument, if 
A ■ U —*■ V is injective, where U and V are open, and K C U is compact, then A -1 , restricted to 
A(K), is continuous. The compactness assumption is one of the ways ill-posedness of many (but 
not all) inverse problems manifests itself. The continuity of A^ 1 us viewed as stability, and such 
kind of results are called conditional stability. One of the central problems in Inverse Problems is 
to estimate the modulus of continuity of A' 1 (restricted appropriately). In other words, we want 
to find a function 4>(t), t — » 0, so that 4> ~ > 0, as t — > 0, and 

||/i-/2|| Si <^(IWi)-^(/2)||b 2 ) 

for all /1, f2 in some neighborhood of a fixed fo, and possibly also, restricted to a compact subset 
K of B\ (hence the stability is conditional). The estimate is of Holder type, if one can choose 
4>(t) = Ct a , < a < 1, and of Lipschitz type, if <p{t) = Ct. We view in this paper inverse problems 
that allow Holder estimates as "stable" (Lipschitz stability occurs rarely); while any other are 
viewed as "unstable". Recall that for the EIT problem, <p(t) = C (log(l/t)) - ^, [x > 0, that tends 
to 0, as t — > 0, much slower than t a . 

The choice of the compact set K also matters. If it is "too small", the conditional stability is 
trivial and not interesting. For example, if K is finite dimensional, and A is smooth enough, one 
always has a Lipshitz stability estimate under the assumption that A has an injective differential, 
see below. If B\ is H S (M) or C k (M) with a compact set (or manifold) M, we are interested whether 
one can choose K to be H 31 (M) or C kl (M) with s\ > s, or k\ > k. A subspace K consisting of 
real analytic functions only, for example, will be considered "too small" . 
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2.1. "Inverse Problems" in R n . We start with the trivial case when B\ = R n , B 2 = R m - Then 
the problem is reduced to solving ^ with A : R n — > R m . Assume that A G C 2 {A G C 1 actually 
suffices, see Theorem Q] below) . Let A XQ be the differential of A at x = Xq. Then 



(4) A(x) = A{xq) + A XQ (x - xq) + R XQ (x) with \R XQ (x) | < C XQ \x - x 



2 



for x near xq. Assume now that A XQ is injective (which can happen only when m > n). This 
immediately implies the estimate 

(5) \h\ < C \A Xo h\, Vh G R n . 

The easiest way to prove © is to note that one can choose 1/Co to be the minimum of \Ah\ for 
h on the (compact) unit sphere. Then the local uniqueness question about the number of possible 
solutions of A{x) = A{xq) near x = xq can be answered as follows. We have 

\x - xq\/Cq < \A X0 (x -x )\< \ A(x) - A{x )\ + \Rx (x)\ 

< \A(x) - A(x )\ + C X0 \x - x \ 2 . 

Therefore, if \x — xq\ < (2CqC Xo )~ 1 , one gets 

(6) \x - x \ <2C \A(x) - A(x )\. 

We therefore have local uniqueness and a Lipshitz stability estimate. Moreover, one can choose C xo 
to be locally independent of xq, if A G C 2 , and the same applies to Co- Thus one gets that there 
exists a neighborhood U of xq in R n , so that for any x\ G U, X2 G U, the equality A(x±) = A(x2) 
implies x\ = x^- Moreover, similarly to ([6]), one has 

(7) \x x -x 2 \< 2C \A(x 1 ) - A(x 2 )\, x x Elf, x 2 G U. 

Note that B 2 can be infinite dimensional here, and the arguments do not change. In particular, 
if A G C l4 \ [i > 0, and A xo is injective, then one always has local uniqueness and the stability 
estimate ([7|) of Lipshitz type. If B\ is infinite dimensional, then A\k has that property for any 
finite dimensional K. 



2.2. A well posed inverse problem. The local injectivity theorem in Banach spaces. We 

return to the case where B\ )2 are finitely dimensional Banach spaces. Consider now a truly mul- 
tidimensional version of our abstract inverse problem under assumptions that make it well-posed, 
and therefore, not typical. The main assumption is that A xo : B\ — > A X0 {B\) is an isomorphism, 
i.e., that the differential A xo is not only injective but also has a closed range. Then we have the 
following theorem that is closely related to the inverse function theorem and the implicit function 
theorems in Banach spaces (A xo is assumed to be an isomorphism from B\ to B 2 then). 

Theorem 1 (Local Injectivity Theorem, PQ). Let A : U C B\ — > B 2 be C 1 , with U 3 /o open, let 
Af be injective and let it have a closed range. Then there exists a neighborhood V C U of fo, on 
which A is injective. Moreover, the inverse A' 1 : A(V) — > U is Lipshitz continuous. 

This theorem shows that we have the same conclusions as in the finitely dimensional case, with 
a Lipshitz stability of non-conditional type. 

We refer to pQ] for the definition of the differential Af of A at a fixed /o and for a definition of 
(a Frechet) differentiable map. In particular, if the Gateaux derivative (the directional derivative 
in all directions) exists in some open set and it is continuous there, then A is Frechet differentiable 
as well and the two derivatives coincide. Then A is said to be C l . 
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2.3. An abstract ill-posed inverse problem. We turn out attention now to an abstract inverse 
problem that may not be well posed. In order to study the inverse problem ([I]) by linearization, 
we assume that the differential Af of A exists at / = /o with some /o- Actually, we will assume 
something slightly stronger — a quadratic estimate on the remainder, that roughly speaking means 
that A is C 1 ' 1 . Namely, we assume that 

(8) A(f) = A(f ) + A f0 (f - /o) + R fo (f), with \\R fo (f)\\B 2 < C fo \\f - f \\ 2 Bl , 

for / in some neighborhood in Si of some /o € B\. 
It may happen that we have the weaker estimate 

(9) II%(/)IIb 2 <C/oII/-/oII1 1 

for / in some neighborhood in B\ of some /o G Si, where Bi C Si, B2 3 S2. Here and below, 
B\ C Si, B2 Z> S2 mean also that 

( 10 ) II • IIbi < c\\ ■ n^, || ■ \\g 2 < C\\ ■ || B2 . 

In such case, we can simply replace S2 by S2; and Si by Si. In the boundary rigidity problem 
[12] . actually Si = C X (M) D Si = C 2 (M), i.e., ([8]) is even stronger. Keeping Si, S2 and assuming 
appropriate interpolation estimates would improve a bit the exponent /ii//2 in the Holder exponent 
in Theorem [2] below but for the sake of simplicity, we will not pursue this. In other words, we have 
some freedom how to chose our spaces Si^, and we choose them in a way that makes ([8]) true. This 
choice might not be the optimal for the analysis of Af however. 

Assume now that Af is injective. This does not imply automatically a stability estimate of the 
type ([5]), of course. Assume in addition, that such an estimate holds in some norms, i.e., 

(11) \\h\\ B{ <c\\A fo h\y 2 , v/iGSi, 

with some Banach spaces B' 2 C S2, B[ D B%. In applications, we want not only such B[, B' 2 to exist 
(they always do, one can set B[ = Si, ||^4y /i||^ = \\h\\si on Ran 4^) but we also want they to be 
"reasonable" spaces independent of /o, typically some H s or C k spaces. Then, if B\ = B[, B2 = B' 2 
one can apply Theorem Q] to prove local injectivity and a Lipshitz stability estimate of the type ([7]). 
The inverse problem is well-posed then. 

If the pairs of spaces in ([8]) and (jlip cannot be chosen to be the same, and often that is the 
case, one can still prove local uniqueness and a stability estimate but of a Holder and conditional 
type, if certain interpolation estimates are satisfied, as we show below. Then the inverse problem 
is ill-posed but we think of it as mildly ill-posed. 

Before formulating this as a formal statement, we will give an example where (jllj) does not hold 
(in any Sobolev spaces), and there is no local uniqueness. 

Example 1. Let I 2 be the Hilbert space of sequences x = {x/ c }^ =1 with norm ||x|| 2 = ^|xfc| 2 . The 
Sobolev spaces h s are defined through the norms ||a;||^ s = fc 2s |xfc| 2 . We think of the components 
of x as the Fourier coefficients of a 27r-periodic function, which explains why we call h s a Sobolev 
space. Define the non-linear map A by 

(12) A(x) = Ex - (x,a)x, 
where 

a = {l/k}, E = di&g{e- k }. 

Clearly, .4(0) = 0, the linearization of A near x = is Aq = E, and the latter is an injective 
map. Moreover, there is an inverse E~ l with a dense domain, but E~ l is unbounded as an 
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operator from any Sobolev space to any other one, thus (jlip does not hold in such spaces. Since 
[*4(x)] fc = e' k x k - x k x m /m, we get that 

x {k) = (0,...ke- k ,0,...), 

where the non-zero entry is the k-th one, is a solution to A{x) = 0. Therefore, in any neighborhood 
of in I 2 there are infinitely many solutions of A{x) = despite the fact that the linearization 
Aq is injective. Moreover, for any s, in any neighborhood of in h s , there are still infinitely 
many solutions. Therefore, there is no local uniqueness in this case, neither weak nor strong, in 
any Sobolev space. One can still get local uniqueness for <C 1 by choosing an h°° type of 

norm || • ||* with an exponential weight, namely = ||22 _1 a;||. This however translates into a 
neighborhood of the origin of 27r-periodic functions that involves certain real analytic functions 
only. We think of such topology as "unreasonably restrictive". 

Based on that example, one can consider the following map in L 2 (R n ) 

A(f) = <f>*f-(f,a)f. 

Choose = Yl X{/c-i<|£|<fc}(£) e_fc '^! where xk stands for the characteristic function of K. Also, 
fix a G L 2 with |a(£)| > (1 + |£|)~ m f° r some m. Then / = solves A(f) = 0, and A has an injective 
differential at / = 0. Also, if / is supported in {k — 1 < |£| < k}, and / solves (f,a) = (27r) n e~ fc l^, 
such an / would also solve A(f) = 0. This gives us a sequence of such solutions, converging 
exponentially fast to in any H s . 

Our main abstract theorem for linearizing inverse problems with a "stable" linearization is the 
following. 



Theorem 2. Let A be as above, and assume |2jj, f77]) with B\ C B[, B' 2 C B2 as above. Assume 
also that there exist Banach spaces B'£ C B' 2 , B'[ C B\ so that Af continuously maps B'[ into B' 2 ' 
and the following interpolation estimates hold 

(13) \\u\\ B r 2 <C|M|g 2 2 H|^ 2 , \\h\\ Bl < CH/ill^ll/ill^ 1 Hi, 112 £ (0,1], MlM2>l/2- 

(a) For any K > there exists e > 0, so that for any f with 

(14) ||/ - /o|| Bl < 6, \\f\\B<><K, 
one has the conditional stability estimate 

(15) 11/ -/oik <C{K)\\A{f)-A{h)r B f\ C{K) = CK 2 ~^\ 

In particular, there is a weak local uniqueness near /o, i.e., if A(f) = A(fo), then f = fa. 

(b) Assume in addition that there is a Banach space K C B'{ so that ill]) holds for fo replaced 
with f close enough to /o in K,, and Af : B'{ — > B 2 ' is uniformly bounded for such f . Then there 
exists e > 0, so that for any f\, fi with 

(16) ||/i-/o|k<e, H/2 -/oik < e, 
one has the conditional stability estimate 

(17) n/i-/ 2 ik<ciwi)-^(/2)iigr- 

In particular, there is a strong local uniqueness near fo, i.e., if A(f±) = A(f2), then f\ = f 2 . 



6 



P. STEFANOV AND G. UHLMANN 



Remark 1. Typical choices of the Banach spaces above are C k or H s spaces. Interpolation estimates 
for those norms are well known. In R", for example, one has 

(18) \\u\\h s < C||it||# si ||u||^ 2 S2 , s = si«i + s 2 a 2 , at\ + = 1, «i > 0, a 2 > 0. 

This can be easily verified by using the Fourier transform. Similar estimates hold for Sobolev spaces 
and the C k spaces in domains/manifolds, see [18]. For example, if M is a compact manifold with 
boundary, then 

(19) IMIo(M') < C 'll n llc 1 ' i i(A/)ll u llc fc 2(M)' k = k 1 a 1 + k 2 a 2 , ax + a 2 = 1, ai > 0, a 2 > 0. 

Interpolation estimates that include an H s and a C k norm can be obtained by using Sobolev 
embedding theorems. 

Remark 2. Assume, for example, that all spaces above are Sobolev spaces. Then one can make 
the exponent [i\\x 2 in (|15p . (|17p arbitrary close to 1, thus making the stability estimate "almost 
Lipshitz". The price to pay for that is to assume that / is a priori bounded in H s with s ^> 1, 
see (I18p . Alternatively, one may try to choose B'{, K, close to £>i, which is the natural limit but 
that will decrease \x\. Since \i\ > l/(2/x 2 ) > 1/2, this puts a limit on how close B'{, K, can be to 
B\. If B\ = H s , B[ = H s , for example, s' < s, then we will have at least B'{ = H s with s" so 
(s' + s")/2 > s. Therefore, we cannot take s" > s arbitrary close to s unless s' = s. The choice of 
s" in this case is further restricted by /i 2 - 

Remark 3. The second inequality in (I14p is a typical compactness condition, when we work with 
Sobolev spaces or C k spaces in a bounded domain or on a compact manifold. It may seem strange 
that there is no such explicit condition in (b). It is actually there, since (|16p implies such a condition 
with B'2 replaces by /C. Actually, using interpolation estimates, one could write (I14p in the form 
(|16p . and vice -versa, with slight change of the Banach spaces (i.e., small change of k in C k , etc.). 
We formulated (I16p in a form different than (|14p in order to avoid introducing yet another Banach 
space. What makes (|15f> and (|17|) conditional estimates is that in the conditions (|14p . respectively 
(|16p . there is at least one norm stronger that the norm in the l.h.s. of (|15p . respectively (|17p . Note 
also that (I15p . (|17p can be formulated with different choices of the norms. 

Proof of Theorem^ We start with (a). By % and CED, 

P/o(/ - /o)b 2 < P(/) - A(f )\\ B2 + C\\f - M\ 2 Bl 
By (fTTj) and the Holder inequality (a + 6) M < a M + 6^ for a > 0, b > 0, we get 

U h {f - fo)\y 2 < CK 1 -^ (\\A(f) - A(fo)\r B 2 2 + 11/ - /oils?) • 

The constant CK 1 ' 112 comes from estimating the term ||^4j (/ — fo)\\ 1 ^, >12 that is bounded by the 

assumption on Aj and by the second inequality in CO}. By (fTTj) . and the second interpolation 
inequality in (fT3T) . 

11/ - /oik < CK 2 -^ (\\A{f) ~ -4(/o)||£f 2 + 11/ - /o||^ lM2 ) • 

Since 2/ii/x 2 > 1, for e <C 1 we get (fT5l) . Note that the condition on e has the form e < 
(C-fT 2 ^ 1- ^ 2 ) -1 , where C depends on .A only. This proves (a). 

To prove (b) we note that the same proof can be applied under the assumptions of (b). Indeed, 
we start with 

11% (f 2 - A) < \\A(f 2 ) - A(h)\\ B2 + C\\f 2 - f x \\l v 
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Then we proceed as above under the condition ||/i — /o||ac < e i ^ 1 that would guarantee (fTT|) 
with /o replaced by f\. Then (flTl) holds if H/2 — /i||bi < ^2 -C 1. Those two conditions certainly 
hold if e in (|16p is small enough. □ 

Remark 4. The variety of norms in the theorem above is typical. For example, in the boundary 
rigidity problem, we first apply A*j q first (called there I*, and the adjoint operator is taken in 
1? spaces with respect to suitable measures). Then Aj is replaced by A*? o Af , and we made the 
following choices, see Hoi: 

(20) B l = C 2 {M), B[ = L 2 (M), B'[ = H S (M), s » 1, K = C k (M), k » s, 

(21) B 2 = L°°{M 1 ), B' 2 = H 2 {M A ), B' 2 " = H S -\M X ), 

where M] nt 3) M, and H 2 (Mi) is an appropriate Hilbert space such that H 2 {Mi) C H 2 {M{) C 
H l (M\). There are other complications there, coming from the fact that we have to work with 
equivalence classes of isometric metrics. Then /*/ acting on symmetric tensors in H S (M), extended 
as zero to Mi, is not a bounded operator with values in H s ~ 1 (Mi) even though 1*1 is a if DO of 
order —1. The reason is that / G H S (M), extended as zero, belongs to H S (M\) only if the first s 
normal derivatives of / (assuming that s is integer) vanish on dM. This is the reason for introducing 
the space H 2 . A stable recovery of the derivatives of the metric on dM from the boundary distance 
function is used in [15] to deal with this problem. 

3. Stability of linear inverse problems 

Most of the material in this section is based on the theory of elliptic ^DOs and on the theory of 
Fredholm operators. We give proofs for the sake of completeness of the exposition. 

Let A : B\ — > B2 be a bounded linear map between two Banach spaces. Assume in addition, that 
A is injective. We are interested whether one has the stability estimate 

(22) < C\\Af\\ B2 

with some constant C > 0. It is well known, that if A is smoothing, and B\ and B2 are among C k 
or H s with finite k, or s, then there cannot be such a stability estimate. We will quickly review 
those arguments below. 

First, note that (1221) is not equivalent to invertibility, because (the closure of) the range Ran^4 
of A can be much smaller than #2. One the other hand, one can replace B2 by B' 2 = Ran A, and 
then (f22|) is equivalent to the invertibility of A : B\ — > B' 2 , as the next proposition shows. 

Proposition 1. 122^ holds if and only if A is injective and Ran A is closed. 

Proof. Assume that Ran A is closed and A is injective. Then A : B\ — > B' 2 has range that coincides 
with B 2 . Since A is injective, then A -1 exists with domain the whole B 2 . By the open mapping 
theorem, A'" 1 is bounded, hence (|22[) holds. 

Assume (|22p . Then the injectivity and the closedness of Ran A follow easily. □ 

Proving that Ran A is closed is not always straightforward. Proving that it is not (when it is 
not), is usually easier. As an example, suppose that Bi = B 2 = L 2 (9), where is a domain in R" 
with smooth boundary, or Q = R™. If there exists s > so that Af € H S (Q) for any / € L 2 (Q), 
then (|22p cannot be true. Indeed, if we know that A has a trivial cokernel (for example, because 
it is injective and self-adjoint, as the operator N w below), then Ran A cannot be closed. In the 
general case, assume that Ran A is closed. Then the latter is a Banach space as well, and A is 
invertible on it. But A is compact, so we get that the identity on L 2 (ft) is compact, that is not 
true. 
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On the other hand, ([22]) may still hold with B 2 = H S (Q), if Ran ,4 is closed in H S (Q). Take for 
example A = (1 - A)" 1 on L 2 (R n ), and s = 2. Then one can ask the more general question: is 
there a stability estimate of the type 

(23) ||/||x2 < C\\Af\\„. 

for some s > 0? If so, it will be a weaker substitute for (|22|) . if the latter does not hold but it will 
still allow us to apply Theorem [2l For (j23[) to be true, A must map L 2 into H s , therefore, this is 
the same question as before, with B 2 = H s . 

Example 2 (the weighted X-ray transform). In R n , let I w be the weighted X-ray transform 

(24) I w f(x,9)= [ w(x,9)f(x + t6)dt, x E R n , 6 E S n ~\ 

where w is a continuous function. Let us restrict I w to functions supported in Cl, where Q C R™ 
is a strictly convex domain with a smooth boundary. One can assume then that (x,8) E d-SQ, 
which means that x E dQ, and 6 E S™^ 1 points into We equip d-S^l with the measure 
dE = \9 ■ u(x)\ dS x d6, where dS x is the measure on dft induced by the Euclidean measure in R n , 
and dO is the measure on S n ~ l . Then one can show that I w : L 2 (S7) — > L 2 (d^SQ, dE) [8], see also 
the remark following (j27[) below. Let 

(25) N w = T W I W : L 2 (n) ^ L 2 (Q) 

be the "normal operator". We can think of the lines through 0, as a manifold, and d-SQ then is 
a global chart for them (singular at the lines tangent to dfl). The measure dE remains invariant 
if we replace <9f2 by another strictly convex surface encompassing S7. Thus one can expect that N 
will not depend on dfl. Indeed, on can easily get that (see [TT] ) 

(26) N w f(x) = c n / iz—r-dy 

J \x-y\ n 1 

with 

_( x — y \ ( x — y\ _( x — y\ ( x — y \ 

(27) W{x,y) = w[ x,— f \w[ y,—. , +w[ x,- , }w[ y,- : . 

V \ x -y\J V \x-y\J V \x-y\J \ \x-y\J 

Now we think of N w as the operator defined by the formula above, with w(x, 6) extended for all x. 
It is easy to see that N w maps L 2 Q (R™) into L 1 2 oc (R"). In particular, N w is bounded on L 2 (Q). 
That shows that I : L 2 (Q) -> L 2 (dSVL, dE) is bounded. 

Let Q,i D 17 be another domain with the properties of O. Then N : L 2 (ft) — > L 2 (Qi) is also 
bounded, and I is injective if and only if N : L 2 (Q) — ► L 2 {Q,\) is injective, where one can replace 
fli by ft as well, see Lemma [1] below. 

If w E C°°(f2), then is a ^I/DO of order —1, as we will see below. Thus ./V maps L 2 (tt) into 
H 1 ^) and a stability estimate for N of the kind ([23]), with B 1 = B 2 = L 2 ({1), or B 2 = L 2 ^), is 
not possible. If w E C 1 ^) only, then one can still show that N maps L 2 (Q) into i? 1 (0) using the 
theory of operators with singular kernels, see [5] for details. We then change B 2 and replace it with 
H l (M\). The question we ask now is the following: assuming that / (and therefore N) is injective 
for some class of w's, what conditions on w would guarantee the stability estimate 

(28) ||/||x 2 (Q) < C\\N w f\\ Hl{ni) , V/ E L 2 (n). 

We will show below that a sufficient condition for that is that ./V is elliptic, that can be easily 
formulated in terms of w. 
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We return to the general case and the estimate (|22p. In the example above, if N w is elliptic, 
then it is also a Fredholm operator. Recall that F : B\ — > B2 is called Fredholm, if there exist 
bounded operators A : B\ — > B2 and B : B2 — > B\ so that — Id and FB — Id are compact. Then 
Ran F is necessarily closed, and F has a finitely dimensional kernel Ker F and a finitely dimensional 
cokernel Coker A = B2/ Ran A . Equivalently, F is Fredholm if its kernel and cokernel are finitely 
dimensional. 

In case the inverse problem is over-determined, we cannot expect that A would have a finite 
cokernel, as in the example above. For this reason, if B2 is a Hilbert space, we will study whether 
N = A* A is a Fredholm operator, and the consequences of that. In the applied literature, applying 
A* is viewed a "back-projection" that, in inverse boundary value problems, returns our data form 
the boundary back to the domain. Another reason to study A* A instead of A is that, as in the 
example above, A might be an FIO, but A* A is a ^DO, thus easier to study. The latter is not 
alway true for any FIO, so some conditions are needed. In case of integrals along geodesies, this is 
guaranteed by the no-conjugate points assumption, that is certainly true for the Euclidean metric. 
Again the example above shows that we may need to restrict N to a subspace £ of H in order to 
avoid working with \&DOs on manifolds with boundary. Then F\c may not have finite cokernel 
anymore. For this reason, we will only study operators F that are upper semi- Fredholm, i.e., Ker.F 
is finite dimensional, and Ran F is closed. Equivalently, F is upper semi- Fredholm, if there exits a 
bounded left inverse A modulo compact operators (AF — Id is compact). Note that if N = A* A is 
Fredholm, then A itself is upper semi-Fredholm. 

We are not losing generality when replacing A by N = A* A as the next lemma shows. Note that 
also stability estimates for N of the type (f28l) and (|29|) below imply stability estimates for A of the 
type (P2}, P3D as well. 

Lemma 1. Let A : B\ — > B2 he bounded, where B\^ are Hilbert spaces, and let C C B\ be a (closed) 
subspace. The following statements are equivalent. 

(a) A : C — > £>2 is injective; 

(b) A* A : C —* B\ is injective. 

(c) licA* A : C — > C is injective, where lie is the orthogonal projection to C 

Proof. Let f £ C. Then (A*Af,f) = \\Af\\ 2 . This proves the implication (a) (b). If Af = 
with / € C, then A*Af = as well, therefore (b) (a). The equivalence of (c) and (a) follows in 
the same way. □ 

Proposition 2. Let N : B\ — > B2 be upper semi-Fredholm. Assume in addition, that N is injective 
on some (closed) subspace £ cTt. Then: 

(a) There exists C > 0, so that 

(29) ll/llft<CW|| fll , V/e£. 

(b ) Estimate \29\) remains true for any other bounded operator close enough to N with a uniform 
C>0. 

Clearly, one can always take C = (KeriV)^ if B\ is a Hilbert space. The way we formulated 
the proposition above however is intended to underline the fact that if we can prove injectivity 
on some explicitly given subspace that is of interest to us (L 2 (f2) in the example above), then we 
automatically have stability. 

Proof of Proposition^ The operator F : C — > F{£) is bounded, injective, and has a closed range. 
In particular, it is also upper semi-Fredholm. By Proposition [H (|29p holds. This proves (b). 



10 



P. STEFANOV AND G. UHLMANN 



To prove (b), it is enough to write 

(30) WfhJC < \\N q J\\ B2 < \\(N q - N qo )f\\ B2 + ||JV,/|| ft . 

□ 

Remark 5. In the tensor tomography problems, one has to use a finer version of Proposition [2](b). 
Then N = N g depends continuously on the metric g, and N g : ff 1 (M) — ► H 2 (Mi), where Mi §> M 
and tensors in M are alway extended as outside M. The space H 2 {Mi) has the property that 
H 2 (Mi) C H 2 {Mi) C H x {Mi). On the other hand, H 2 (M{) is "too large", and R&nN g is not 
closed there. It turns out, that one can construct a left parametrix Q g so that Q g N g = Id + K g , 
where K g is a smoothing operator. However, Q g : H 2 (M\) — ► H l (M) is not bounded; we can show 
only that Q g : # 2 (Mi) -> L 2 (M) is bounded. On the other hand, Q g N g is bounded on H l (M), of 
course. This requires some modification of the arguments above, and the resulting estimate is 

\\f s \\ LHM) < C\\N g f\\ kHMy V/ € H\M) 

with C > locally uniform in g. Here / s is the orthogonal projection of / onto a certain space, 
called solenoidal tensors, where N g is injective by assumption. Note that (f30|) does not apply with 
B x = L 2 {M) and B 2 = H 2 {Mx) because N g is not bounded as an operator between those two 
spaces. On the other hand, one cannot replace B\ by H (M) because N g (H (M)) is not closed in 
H 2 (M). For more details, and for the needed modifications in the arguments, we refer to I12j. 

The following proposition is well known, see e.g., |17j . We define Sobolev spaces on a compact 
manifold by fixing a finite atlas. We assume in the proposition below that N is of order — m < 
in view of applications to tomography problems. The order —m does not need to be non-positive, 
and one can study iV as an operator mapping H s into H s+m for s not necessarily by the usual 
technique — applying powers of 1 — A to the left and right. 

Proposition 3. Let N € *$>~ m (M), m > be an elliptic ^>DO on a compact manifold M . 

(a) Then for any I > 0, one has 

(31) \\f\\l?(M) < C, (||A7||^(M) + ll/lk-HM)) • 

(b) Assume in addition, that N is injective on some (closed) subspace C C L 2 (M). Then there 
exists C > 0, so that 

(32) ll/lb(M) < C\\Nf\\ Hm(M) , V/ G C. 

Proof. Part (a) follows directly from the existence of a parametrix. Indeed, there is Q € ^> m so 
that QAf = f + Kf, where K has a smooth kernel. Note that for this proof it is only enough 
to have a parametrix Q of finite order I, i.e., to have K that maps L 2 into H l . Then flS]) follows 
directly. 

Note that the inclusion L 2 — > H~ l is a compact map. Part (b) follows from Proposition [2] above. 
Note also that an alternative way to prove (b) is to notice that it follows directly from (a), by 
[171 Proposition 5.3.1]. □ 

Recall that Vl/ — m (M) is a Frechet space with semi-norms pk(A), k = 1,2,... given by the semi- 
norms (they are actually norms) of its amplitude in the finitely many charts (Uj,Xj) of M. The 
norms pk are given by 

Pfe (a) = max sup max (1 + Id^dja^x, y, 

3 se & f ,€6R»l a l+l/ J l+W<* 

where aj is the full symbol of A in the j-th coordinate chart. 
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Proposition 4. There exists k 3> 1, depending only on dimM with the following property. Let 
A = Aq G ^>~ m satisfy the assumptions of Proposition \2(b). Then 129\) holds for any A G fy~ m 
close enough to Aq in the norm pf. with a constant C independent of A. 

Proof. Note that <J/ _m equipped with the norm p^, k > 2n + m + 1, is a subspace of the space 
of linear bounded operators mapping L 2 (M) into H m (M), see the appendix in [10]. Now the 
proposition follows from Proposition Wijo). □ 

Remark 6. Proposition [3] can be applied to compact manifolds M with boundary; then we require 
the symbol to be defined and to satisfy the symbol estimates for x in the closed M. One can then 
extend it as an elliptic symbol to some extension of M without boundary, and treat L 2 (M) as C. 
The same remark applies to Proposition |4] because on can choose the extension of the symbol in a 
way that is continuous in any fixed pk norm. 

Finally, we formulate a direct corollary of these statements. 

Corollary 1. Let A : L 2 (M) — > TL, be a bounded linear operator, where M is a compact manifold 
without boundary, and TL is a Hilbert space. Assume that N = A* A is an elliptic ^DO of order 
—m < 0. Assume also that A is injective on a (closed) subspace C of L 2 {M). 
(a) Then N : C — » C is onto, and 



(b ) Estimate A33\) remains true with an uniform C > if N is replaced by any other operator 
close enough to N in the operator norm L 2 (M) — > TL, and in particular for all ^DOs in a small 
enough neighborhood of A in the norm pk, k S> 1 . 

Proof. The only statement left to prove is that N : C — > C is onto. This follows from the fact that 
N is injective and self-adjoint, hence it has a trivial cokernel. Next, since N is Fredholm, its range 
is closed. □ 

In particular, by Lemma [U Corollary [TJ^b) implies that A remains injective under a small per- 
turbation (which also follows from the fact that A\c is upper semi- Fredholm) . Also, (|33p implies 
an estimate of the kind ||/|| < C||j4/||* with a suitable norm || • ||*. 

Note that Corollary [T] can be applied to manifolds with boundary as in the example below. 

Example 3. We return to Example [21 Let £l\ D f2 as before. Assume that N w = L^I W is an elliptic 
^DO (of order —1) in Q\. To apply Corollary [1] directly, one can embed f2i into a compact manifold 
and extend N w as an elliptic operator there in a way independent of w. Alternatively, on can work 
with the parametrix in f2i and eventually restrict to £1. By (12TH . the principal symbol of N2 is 
given by (see [HI [I2j [5] derivation in even more general cases) 



where 5 is the Dirac delta function, and clearly, a p (N w ) is homogeneous in £ of order —1. The 
singularity at £ = can be cut-off. Since w is at least continuous, we see that the ellipticity 
assumption is equivalent to the following: 



(33) 



f\\mM) < c\\Nf\\ Hm(M) , v/ g c. 




(34) 



v(x, c) g n x S' 



tn—l 



39 G S' 



m— 1 



9 _L C; so that w(x,9) / 0. 



In invariant terms, (x, £) G S*fl, while (x, 9) G Sfl. Another way to express the same is to say that 
(35) iV*(supp u>) D T*0, 
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where supp w is the open set where a) / 0. Under this condition, using the continuity of w in 
SQ = f2 x S n , one can get the same with Q replaced by fli. By Proposition [31 we then get 

(36) \\f\\ L 2(u) < C\\N w f\\ m(Ql) , 

under the assumption that I w is injective on L 2 (0). Since N w is elliptic, it is enough to assume 
only that I w is injective on Cq°(Q,). Without that assumption, one has fj31[) : and ()36p on (KeTl w )~ L . 

To analyze the continuous dependence on w, let us note first that given k > 0, one can construct 
a bounded extension operator extending functions in C k {Cl) to C k (Cli). With the aid of this 
operator, one can map a neighborhood of some wq £ C k (Cl) to some neighborhood of an extended 
wq in C k (Qi). Then we can apply Corollary [1] to get the following. 

Theorem 3. Let w € C l {pt x S n ~ l ) satisfy ( f5^| ). Assume that I w is injective on C^°(r2). Then 
N w : L 2 (fi) — > L 2 (Q) is onto, and there exists a constant C > so that 

(37) WfWvm < c\\N w f\\ L2{ni} , v/ g L 2 (n). 

Moreover, (31) remains true under small C 1 (fi) perturbations ofw. 

Without the injectivity assumption, (|37p holds on (Ker I w \cg° (Q))~ L ■ Theorem [3J is a weaker 
version of the results in [5] and is presented here just to illustrate the method. Note that w € C k 
with k = 1 is not enough to guarantee that operator convergence in the Pk{N) norm implies 
convergence in the norm H-iV^^w^ffJi)) viewing N = N w as a ^DO. On the other hand, N w is 
also an operator with singular kernel, see (|27|) . and w G C 1 is enough for I? — ► ^T 1 continuity. We 
refer to [5] for details. 

We show finally that for an injective ^DO, ellipticity is a necessary and sufficient condition for 
the estimate (j33j) in Corollary CD to hold. 

Proposition 5. 

(a) Let N be a ^>DO of order —m < on a compact manifold M. Then A32\) holds for any 
f £ L 2 (M) if and only if N is injective and elliptic. 

(b) Let M be a compact manifold with boundary and let M\ be another compact manifold so that 
Mj nt D M. Assume that N is a ^?DO of order —m < defined on some neighborhood of M. Then 

(38) \\f\\LHM)<C\\Nf\\ Hm{Ml) , V/€L 2 (M) 

if and only if N is injective on L 2 (M) and elliptic on tt^ 1 (M); more precisely, 

inf hminf|erbo(x,e)l>l/C, 

xeM too 

where po is the principal symbol of N . 

Proof. To prove (a), we will work locally near a fixed point on M. In fixed local coordinates near 
xq, let p(x, £) be the symbol of N, multiplied by a standard cut-off supported near some xq. We 
will study p(x, D) in R n first. Fix ^ and apply N to the normalized "coherent state" 

f x ( X ) = ( A / 7r )n/4 e iAx. € o-A|x~xo| 2 /2_ 

Then 

/ A (A£) = {\/n) n l\2ir/\) n l 2 e- iX ^y^-^-^\ 2 /z. 

Then we get 

\p(x,D)f x ](x) = e iAl, °- & (A/vr) n / 4 (27rA)- n/2 f e iXix - x ° H - x ^-^ 2/2 p(x, \£) d\£. 
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For any C > 0, the contribution to the integral above coming from integrating over |£ — £q| > 1/C 
is 0(A ), ViV > 0. To estimate the effect of replacing p(x, At/) by p(x, A£o) in the integral above, 
notice first that p(x, A£) — p(x, A£o) = C( x : £> Co ; A) • (£ — Co) with £ = 0(1), therefore 



/ 



(39) / |e-Co|e- A| ^ ()| /2 dC = 0(A"") 

l<c 



(pass to polar coordinates rtu = £ — £q)- Thus we get 
|p(s,D)/ A ](a;) 

= e iAx -5o( A / 2vr )"/2( A / 7r )n/4 



e iA(x-z„K-A|{-Co| fipfa A ^ Q ) d ^ + ( A -n) 
e iAx -Co( A / 2vr )™/2( A / vr )n/4 A ^ )( 2 ^/ A )n/2 e iA(a-x )-€o-A|x- a; o| 2 /2 + 0(A~ n ) 

p{x,\io)h{x)+0{\- n ' A ). 



Assume now that (|32p holds. It is enough to consider the case m = 0. Then, choosing / = fx, we 
get 

(40) 1/C< ||p(-,A£ )/A(-)ll+O(A- n/4 ). 

The function </>a(£ ~~ x o) = I/a( :c )| 2 = (A/7r) n / 2 e _A ' :c_:ro ' 2 is normalized so that it has L 1 norm equal 
to 1. Therefore, 

(41) ||p(.,A£o)/a(-)H 2 = J Mx-x )\p(x,XCo)\ 2 dx=\p(x ,X^)\ 2 + 0(X- n ) 

by the argument already used above, see (|39|) . Relations (|40p and (|4ip show that p is elliptic at 
(xo,£o)- More precisely, we get 

(42) ip(x,oi>i/c-Ci[er ,l/4 ) V(i,()erM, 

where C > is the constant C in (|32p . and C\ > depends on j? only. 

To relate N and p(z, -D), note that in a fixed coordinates system near some Xq, R = N — p(x, D) 
is smoothing when acting on distributions supported near xq. Consider then (N — p(x,D))xf\, 
where x ls a standard cut-off near xq. Then 

p(x, D)f x = P (x, D)(l - X )/a +p(x, D) x fx = p(x, D)(l - x )fx + (N — R) x f\. 

Now, (1 — x)f\ = 0(e~ x / c ), Nxfx = 0(A~°°) as can be seen by integrating by parts and using the 
smoothness of the kernel of R. Therefore, Nxfx = p(x,D)fx + 0(A~°°), and this completes the 
proof of (a). 

To prove (b), we proceed in the same way in a neighborhood of any xq € M mt but we apply 
p(x,D) to xfx instead to fx, where x is a standard cut-off near xo supported in M. As xo gets 
closer to dM, the derivatives of x w iU have larger norms. This will affect the constant C\ in (|42p 
that nay not be uniform in x. The conclusion (b) is still true, however. Note that if N has a 
polyhomogeneous symbol, then the sub-principal symbol must be uniformly bounded for x 6 M, 
because we assume that N is a ^DO in the larger M\. Then C\ is actually uniformly bounded for 
x € M. □ 



14 



P. STEFANOV AND G. UHLMANN 



3.1. Non-sharp linear stability estimates. In view of the analysis in Section [21 we are also 
interested whether an estimate of the type (jlip . weaker than (j33[) can still be proven when N is 
not an elliptic ^DO. In other words, what conditions would imply (|33|) in different H s spaces? 
Ellipticity of N is not necessary anymore, and for certain classes of hypoelliptic operators, one can 
still have (|33p with a loss of one derivative, for example. We are not going to review those classes 
of operators. Without a proof, we will only mention that if the full symbol of N vanishes on some 
open conic set, (|33p cannot hold in any Sobolev spaces. This is well known and used in tomography 
to determine subsets of lines/curves on which the X-ray transform can or cannot determine the 
function we integrate in a stable way. In particular, in Example [21 if T* \ N*(suppw) contains an 
open conic set, see ([35]) . then no H S1 — > H S2 stability estimate of the type (f36|) is possible. 

3.2. Injectivity of the differential Af . In this section so far we were trying to find sufficient 
conditions that would guarantee a stability estimate under the assumption that the differential A f Q 
is known to be continuous. It is not the scope of this work to study in detail the typical approaches 
to prove that the differential is continuous. We will only briefly mention some of them: smallness 
of the coefficients is often used when the the constant coefficient case is easy to study directly. 
Another technique is the use of the Analytic Fredholm Theorem [7] that, when can be applied, 
gives injectivity of Af for generic /'s. When the unknown coefficient is in the principal symbol, as 
in the boundary rigidity question and related hyperbolic inverse problems |13j . the calculus of the 
analytic ^DOs can be applied. 

4. An example: The inverse back-scattering problem for the acoustic wave 

equation 

Consider the acoustic wave equation 

(43) (<9 t 2 - c 2 (x)A) u = 0, t G R, x G R 3 
with a variable speed c(x) > that equals 1 for large x, i.e., 

(44) c(x) = 1 for \x\ > p 

with some p > 0. Let S(s,u>,9), where (s,u>, 9) G R x S 2 x S 2 ), be the scattering kernel associated 
with c, see below. In [10] we showed that if c is close enough to 1 in the C 9 norm, then the 
back-scattering kernel S(s,—9,9) determines c uniquely. 

We recall some facts about the time-dependent scattering theory for (|43[) . see |10j . For simplicity, 
we work in space dimension n = 3. Given 9 G S 2 , we denote by u(t,x,9) the distorted plane wave 
defined as the solution of (l43j) satisfying 

(45) n| t « = 5{t-x- 9). 

Then we set u sc = u — 5(t — x ■ 9). In the Lax-Phillips scattering theory [6], see also The 
asymptotic wave profile of u is defined by 

ul c (s,u,9) = lim (t + s)d t u BC (t, (t + s)u>, 9), 

t^oo 

where s G R, u> G S 2 . The scattering kernel S(s,uj,9) is given by 

S(s,uj,9) = -—uI c (s,uj,9). 

Then S(s' — s,ll>',u>) is the Schwartz kernel of 1Z(S — Id)7?. _1 , where S is the scattering operator, 
and 1Z is the Lax-Phillips translation representation. 
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One can formulate this problem with stationary data in mind. The scattering amplitude a, that 
is also defined through the stationary scattering theory, is essentially just a Fourier transform (a 
stable and invertible operation) of S in the s variable, more precisely, 

(46) — a(A,w,0) = [ e- isX S(s,LU,6)ds, 

2ir J 

see [ID]. 

Set ijj = —9. Then 

(47) supp,S(-,-0,0) C (-oo,2p]. 

Let Pq be the plane x ■ 9 = —p, and let "dist" be the distance function in the metric c~ 2 dx 2 . Fix 
T so that 

(48) T > 2max{dist(P e ,x); x G B(0,p), 9 G S 2 } - 2p, 

where B(0, p) is the ball with center and radius p, see also (|59p . We use the information contained 
in S for s > —T only. We will denote by h the Heaviside function (the characteristic function of 
R+). Our data now is 

(49) S T (s, -0, 9) := h{s + T)S{s, -9, 0). 

Instead of working with St, we will work with its stationary analog defined by 

(50) a T {\,-9,9) := ^ e- lsX S T {s,Lo,9)ds, A > 0. 

The extra factor A -2 there, compared to (|46p . is there for convenience. Clearly, \ 2 a,T = St * XT, 
where the hat indicates Fourier transform w.r.t. s. Getting a stability estimate in terms of ay 
instead of a is even stronger because it uses less information. Since ay is analytic in A, we can 
throw away any finite interval without losing information. Our next theorem shows that we are 
not losing stability, either. We fix Ao > and restrict A to the interval A > Ao. 

Theorem 4. There exist e > 0, k > 2, p G (0, 1) with the following property. For any two c, c 
satisfying (44\ ) an d 

(51) \\c — l||c* + ||c — l||c* < £j 
and for any Ao > 0, we have 

(52) \\c - c\\ L oo < C [ sup (l + A)|(a T -a T )(A,-6>, 

\A>A ,6»e5 2 

where T is fixed so that it satisfies |^<g[ ) w.r.t. both c and c. 

Proof. We start with standard geometric optics arguments, see also |10l Prop. 3.1]. Let hj(t) = t 3 / j\ 
for t > 0, and hit) = for t < 0; then hj = hj—i, j > 0, with the convention = h, h—i = 5. 

Proposition 6. If c G C°° , and eCl, then ViV, 

(53) n= aj(x, 9)hj(t — 4>(x, 9)) + rjy(t, x, 0) 

/or |x| < where <f> solves the eikonal equation 

(54) c 2 (x)|V</>| 2 = 1, 0U. e <_ p = x • 0; 
a,- sofce t/ie corresponding transport equations; and r^ G C N . 
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Proof. The construction of otj is standard, we will focus our attention on rjy. The latter solves 

(d 2 - c 2 (x)A) r N = (c 2 Aa N )h N (t - 4>(x,6)), r N \ t< - p = 0. 

By standard energy estimates, rjy G H^ 1 , see also [ID] . By Sobolev embedding theorems, we get 
rjy £ C^" 1 . To prove that € C , we apply this to rjv+i; then ayvH-i^JV+lC^ — ^>(x, 0)) + rjy + i € 



Note that p + 0(ac, 0) = dist(Pg,x). The expansion ([53]) is valid in B(0, p) if the eikonal equation 
has a smooth solution there. The latter is guaranteed by the following geometrical condition. Let 
Izfiif) be the geodesic in the metric c~ 2 dx 2 issued from z in the direction of 9. Then we want 
z G Pg, and t to be coordinates in B(0,p). The is certainly true if c is close enough to a constant 
m C 2 that we assume. The leading coefficient a_i solves a homogeneous transport equation, and 
in particular, q_i > 0. If c = 1, then a_i = 1, ctj = for j > 0, and <f> = x ■ 9. 

Note that ([53]) remains true if c <E C k , k > 2 but for some iV = iV(fc), and N(k) > 1 if k » 2. 
Moreover, {aj}^_ 1 depend continuously on c G C fc . Therefore, under the assumptions of the 
theorem, for any fixed N, 

(55) ||a_i — 1||(7jv < Ce, ||aj|| c jv < Ce, j = l...N, \\rN\\ C N<Ce, 



if k ^> 1. The norm of ay is taken in B(0,p) x 5 2 , while the norm of r^v is taken in for 
(x,0) G B(0,p) x S 12 , and t in any fixed finite interval (under the assumption e<l). Similarly, 

(56) \\<p - x ■ 9\\ c n < Ce, 

see also [TU] . 

Let c and c be two sound speeds, and denote by S, u the scattering kernel S and the solution u 
related to c. The following formula is proven in [10J: 



(57) (5- S)(s,u,9) = JJ (c~ 2 - c- 2 )u(t,x,9)u(-s - t,x,-to)dtdx. 

The integral above makes sense even though it and u are distributions. 
By the finite speed of propagation for ([15]). see e.g. [3], 

(58) suppu(i, •, 0) C {a;; 0(x, 6) < t} , supp n(—s — t, •, 0) C {x; 0(x, 0) < — s — t} . 
Therefore, for s > —T, on the support of the integrand in (|58p . we have 

t < T-cf)(x,9), -s-t<T-cj)(x,6), 

see ([58]) . This gives an upper bound for t and —s — t for |x| < p. We want that upper bound 
to be so that the singular front of both waves involved in (|46p . with u = —9 there, would pass 
through the whole ball B(0,p). Physically, this is equivalent to the requirement that we make our 
measurements for an interval of time so that the front of the wave u(t, x, 9) would have enough 
time to pass through the whole B(0,p), and then return back a scattered signal. This explains the 
choice of T: we want T to be so that T — 4>{x,9) > (p(x,9), Vx G B(0,p), i.e., 

(59) T > 2max-{>(x,0); M < P, \0\ = 1}- 

This is equivalent to ([38]) . Therefore, we study S restricted to — T < s < 2p (the upper bound 
is not a restriction actually, see (|47[) ). and the length of that interval, assuming for a moment 
equality in (|4"8]). equals twice the length of the longest geodesic issued from some support plane 
Pq = {x ■ 9 = —p}, and staying in 5(0, p). 

For x G B(0, p), s > —T, and t > T + p, we have that — s — t < —p, and then u(—s — t, x, 9) = 0, 
see ([35]) . Similarly, for x € B(0,p), s > —T, and — s — t > T + p, we have that t < —p, and then 
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u(t,x,9) = 0. This shows that in the integral (|57p . the solutions u and u can be replaced by the 
same solutions restricted to t < T + p, and this will not change the l.h.s. when s > —T. Let x De 
a smooth function so that 

(60) X (t) = 1 for -p < t < T + p, X £ C °°(R). 
Set 

(61) u T = x(t)u(t,x,9), 

and we define ut in a similar way. Note that the subscript T in S and u has a different meaning. 
The cut-off for t < —p in ([611) does not change u for x € -6(0, p) because u = there. Then 

(62) (§t ~ S T )(s, -6, 9) = ^ 3 S 3 / ^ (c~ 2 - cT 2 )u T (i, x, 9)u T (-s - t, x, 9) dt dx. 
Take Fourier transform of both sides to get 

(63) (or - or)(A, -0,0) = --^ / (c" 2 (x) - c- 2 (x)) v T (x, 9, X)v T (x, 9, A) dx, 
where 

(64) v T (x,9,X) = J e iXt u T (t,x,9)dt. 

Remark 7. Note that if we replace ut by u on the r.h.s. above, we get the "distorted harmonic wave" 
solution v = e lXx ' e + v sc , where v sc is outgoing, i.e., it satisfies the Sommerfeld radiation conditions. 
In our case, vt is a convolution of v with the Fourier transform of a certain characteristic function. 
That convolution does not affect the high A asymptotic of v, so it is not surprising that below we 
get the expansion familiar from the stationary theory. The advantage that we have here is that 
we do not need to estimate the contribution of large t's, and instead of using estimates of the 
cut-off resolvent for non-trapping systems, we can just use the geometric optics expansion (I53p that 
is more direct, easier to justify for c's of finite smoothness, and easier to perturb near a fixed c. 
Furthermore, we get an even stronger result by discarding information (i 3> 0) that is not needed. 



Linearizing (|63j) . we get the following. 
Proposition 7. Let c, c be two speeds satisfying the assumptions of Theorem^ Then 

(65) (a T -a T )(X,-9,9)= f (c~ 2 (x) - c~ 2 (x)) v T (x, 9, A) dx + R(9, A) 



where 

(66) sup (l + \)\R(9,\)\<C\\c- 2 -c~ 2 \\ 2 Loo . 

ees 2 ,\>o 

Proof. The difference u — u solves 

(a 2 - c 2 A)(n - it) = (c 2 - c 2 )An, (u - n)| K<0 = 0. 

By standard energy estimates, in any compact set, 

\\ut - ih\\ L 2 + \\u - u\\ H i < C\\c 2 - c 2 || loo . 

The constant C depends on c and c and is easy to see that it remains uniformly bounded, if c is 
bounded in C 1 that follows from the hypothesis (I5ip . Then we get by ([61 



(1 + A)||0r -«r||ra <C\\c 2 -c 2 \ 
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Compare ([63]) and ([65]) to get 

sup|(l+A)i?(0,A)l <C\\cT 2 -c- 2 \\ 2 Lao , V0, 

A 

where the L 2 norm is taken in any compact. □ 
We focus our attention now on the linear operator 

(67) Af(X, 9) = J 4(x, 9, X)f(x) dx, f G L 2 (B(0, p)). 
The geometric optics expansion ([53]) yields 

(68) v^(x, 9, X) = e l2X ^ x ^ (a_i(x, 9, X) + ao(x, 9, X) + ...0^(2, — > 00, 
where 

e iXt aj{x,9)x(t)hj{t - (f>(x,9))dt, R N = e iXt x(t)r N (t,x,9) dt. 



Set r/ = 2A#, and extend by homogeneity (of order 1) w.r.t. to 9; then 2X(p(x,9) = <j)(x,rj). The 
leading term is 

a 2 „ 1 = X 2 (0(x,^))Q 2 _ 1 . 

We have that <j)(x,9) is 0(e) close to x ■ rj, see ()56|) . The natural choice for the space containing 
Ran A is therefore L 2 (R"), that we will restrict to \r]\ > 2Ao- In terms of the original variables, 
this space is isomorphic to L 2 ((Ao,oo) x 5 2 ; A 2 dAd#). Clearly, if c G C°°, ([68]) implies that the 
expression in the parentheses on the r.h.s. is a classical elliptic symbol. For c G C k , k 1 the 
symbol satisfies the symbol estimates up to a finite order m only, with m ^> 1, when k 1. To 
prove this, we first show that a.,- is a symbol of order —j (of finite regularity). The properties of r^v 
that we established do not prove that Rn satisfies the symbol estimate \X^ cf^d" e rN\ < CX~ N but 

it is clear that it satisfies \&x^x e r N I < CX~ N . We can replace N by iV + m now and this argument 
is enough to show that A* A is a \&D0 with symbol of finite regularity m. This is summarized in 
the following. 

Proposition 8. Let c G C k (R n ) satisfy T/ien 

vr^/M, M/2) = e iA ^>^a(x,77) 

where 4> and a belong to C N (-B(0, p) x S 2 ) with N = N(k) — > 00, as k — > 00; a(x,rj) is a formal 
elliptic symbol of order satisfying the symbol estimates with all its derivatives up to order N . 

We now consider A as the operator given by (|67|) in the following spaces 

(69) A : B[ := L 2 (B(0,p)) — ► B' 2 := L 2 ((A ,oo) x S 2 ; A 2 dAd#). 
Then 

(70) A*Af{x)= ft f e-' lX2 ^ xfi) -^ y > e)) a{x,9,X)a{y,9,X)f{y)dyX 2 dXd9. 

Js 2 J \>\q J\x\<p 

By a standard argument, define £ = £,(x,y,r]) by 2X<fi(x,9) — <p(y,9) = (x — y) ■ £, where rj = 2X9 
as above. The latter equation certainly has a solution for \x\ < p if c is close enough to a constant 
in C k , k 3> 1, and that solution is close enough to £ = r] in C 1 ' with / 3> 1 when I; > 1. Then we 
make a change of variables from r\ to £ to get that j4*j4 is a \&DO of order in a neighborhood of 
B(0,p), and it is clearly elliptic. 
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To apply the results in section [31 and especially Proposition [3jb), we have to show that A* A 
is injective. Set M = B(0,p\), where p\ > p is close enough to p so that the construction above 
remains valid on M. Then set C = L 2 (B(0, p)). Assume first that c = 1. Then 

A*Af(x)= [ [ [ e~ i2X ^- y > e f( y x)dyX 2 dXde = TT 3 h(\D\-2X )f, 

Js 2 J\>X J\x\<P 

where / is extended as zero outside B(0, p). It is easy to see that h(\D\ — 2Ao) is a Fredholm injective 
operator on L 2 (B(0, p)), see [9] for any Ao > (and it is the identity of Ao = 0). Therefore, the 
conditions of Proposition El^b) are satisfied for c = 1 and we get 

(71) ll/llL 2 (B(0,p))/C < \\ A * A f\\L 2 (B(0,R)) < C\\f\\L 2 (B(0,p))- 

with some R > p. This show in particular that A : B[ — > B' 2 is bounded, depending continuously on 
c G C k , > 1. That property is preserved if we increase p a bit. Therefore, A* : B 2 —> L 2 (B(0, R)) 
is bounded a well, and also depends continuously on c S C k , k 3> 1. Thus 

(72) ll/lb( B (o,p)) < C\\Af\y 2 . 

This is the basic estimate that we need to apply Theorem [2j The spaces B\ and B2 are determined 
by the estimate in Proposition [TJ set B\ = L°°(B(0, p)), and let B2 be the Banach space determined 
by the norm on the l.h.s. of (|66p . We need to show that the interpolation estimates (1131) hold with 
a proper choice of B'[ 2 ■ Since 

ll/llsi = ||/||i°°(B(o,p)) < C e \\f\\ Hn/2 +e, 

Ve > 0, we use (fTHj) to deduce that the second interpolation estimate in (fT3j) holds with = 
H S2 (B(0, p)), where n/2 + e = £2(1 — pi)- Note that we can make pi as close to 1 as we want by 
choosing S2 large enough. Next, we need to choose B' 2 ' so that the first interpolation estimate in 
(USD holds: 

(73) H|^<C||< 2 2 ||u||^ 2 . 

In this case, we can just choose B 2 = £>2, and P2 = 1 because it is easy to see that 

(74) \\u\\ B , 2 <C\\u\\ B2 . 

The condition p\p2 > 1/2 is therefore satisfied when S2 > n + 2e. Then we can choose p = pi to 
be any number in (0,1) so that 1 — p < n/(2s2). The compactness assumption < K, see 

dH]), is then satisfied if ||/||£P2 < K, which is guaranteed by (f5"T|) . if k > S2. 

This completes the proof of Theorem HI □ 
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